Event-by-event background in estimates of the chiral magnetic effect 
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In terms of the parton-hadron-string-dynamics (PHSD) approach - including the retarded elec- 
tromagnetic field - we investigate the role of fluctuations of the correlation function in the azimuthal 
angle ip of charged hadrons that is expected to be a sensitive signal of local strong parity violation. 
For the early time we consider fluctuations in the position of charged spectators resulting in electro- 
magnetic field fluctuations as well as in the position of participant baryons defining the event plane. 
For partonic and hadronic phases in intermediate stages of the interaction we study the possible for- 
mation of excited matter in electric charge dipole and quadrupole form as generated by fluctuations. 
The role of the transverse momentum and local charge conservation laws in the observed azimuthal 
asymmetry is investigated, too. All these above-mentioned effects are incorporated in our analysis 
based on event-by-event PHSD calculations. Furthermore, the azimuthal angular correlations from 
Au+Au collisions observed in the recent STAR measurements within the RHIC Beam-Energy-Scan 
(BES) program are studied. It is shown that the STAR correlation data at the collision energies of 
y/SNN = 7.7 and 11.5 GeV can be reasonably reproduced within the PHSD. At higher energies the 
model fails to describe the ip correlation data resulting in an overestimation of the partonic scalar 
field involved. We conclude that an additional transverse anisotropy fluctuating source is needed 
which with a comparable strength acts on both in- and out-of-plane components. 

PACS numbers: 25.75.-q, 25.75.Ag 



I. INTRODUCTION 

A fundamental property of the non-Abclian gauge the- 
ory is the existence of nontrivial topological configura- 
tions in the QCD vacuum. Spontaneous transitions be- 
tween topologically different states occur with a change 
of the topological quantum number characterizing these 
states and induce anomalous processes like local violation 
of the V and CV symmetry. The interplay of topological 
configurations with (chiral) quarks shows the local imbal- 
ance of chirality. Such a chiral asymmetry when coupled 
to a strong magnetic field induces a current of electric 
charge along the direction of the magnetic field which 
leads to a separation of oppositely charged particles with 
respect to the reaction plane [ll-Q . 

This strong magnetic field can convert topological 
charge fluctuations in the QCD vacuum into a global 
electric charge separation with respect to the reaction 



plane. Thus, as argued in Rcfs. [l|, 0-H|: the topologi- 
cal effects in QCD might be observed in heavy-ion col- 
lisions directly in the presence of very intense external 
electromagnetic fields due to the "Chiral Magnetic Ef- 
fect" (CME) as a manifestation of spontaneous violation 
of the CV symmetry. Indeed, it was shown that electro- 
magnetic fields of the required strength can be created 
in relativistic heavy- ion collisions P, [a, 0] by the charged 
spectators in peripheral collisions. 

The first experimental evidence for the CME, identified 
via the charge asymmetry, was obtained by the STAR 
Collaboration at RHIC at ^/sWn =200 and 62 GeV @- 
[Toj and confirmed qualitatively by the PHENIX Collabo- 
ration [Tl| . Recently, these measurements were extended, 
from one side, below the nominal RHIC energy down to 
^Jsnn =7.7 GeV within the RHIC Bcam-Energy-Scan 
(BES) program fl2l ] and, from the other side, preliminary 
results for the maximal available energy ^/s~nn =2.76 
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TeV were announced from the Largc-Hadron-Collidcr 
(LHC) 00. Though at first sight, some features of 
these data appear to be consistent with an expectation 
from the local parity violation phenomenon, the interpre- 
tation of the observed effect is still under intense discus- 
sion [H-il. 

The fluctuation nature of the CME will give rise to a 
vanishing expectation value of a V-odd observable and 
due to that, as proposed by Voloshin [2]|, the azimuthal 
angle two-particle correlator related to charge asymmetry 
with respect to the reaction plane is measured in exper- 
iments (o-liol. Il2l - fl4| . Accompanying these experiments 
hadronic estimates of the dynamical background in these 
experimental papers including only statistical (hadronic) 
fluctuations do not involve the electromagnetic field at 
all. The electromagnetic field - created in heavy-ion col- 
lisions - was calculated in different dynamical approaches 
inRefs. [3, IE 0, ■ In two of them @,|27j] calculations 

were carried out in comparison with the CME observable. 
However, in all these studies only the mean electromag- 
netic field was presented, being averaged over the whole 
ensemble of colliding nuclei. 

As noted in Ref. [29[ event-by-event fluctuations of the 
electromagnetic field in off-central heavy-ion collisions 
can reach rather high values comparable with the av- 
erage values. The presence of large fluctuations was then 
confirmed in a more elaborated model in Ref. (30j . In 
this study - based on the parton-hadron-string dynam- 
ics (PHSD) kinetic approach - we analyse event-by-event 
fluctuations in the electromagnetic fields as well as the in- 
fluence of these fluctuations on other physics quantities, 
in particular to measurements of the CME. 

The paper is organized as follows: After a short re- 
capitulation of the PHSD approach in Section I we se- 
quentially consider the manifestation of the initial ge- 
ometry fluctuations in spectator protons and participant 
nucleons as well as in the charged quasiparticle geometry 
at some later stage. These effects are relevant for fluc- 
tuations in the electromagnetic fields, the event plane 
orientation and the possible formation of a fluctuating 
electric charge dipolc/quadrupole transient subsystem, 
respectively. Conservation of the transverse momentum 
and local charge is analyzed as an alternative explana- 
tion of the observed azimuthal asymmetry. In Section 
III we discuss the role and importance of these effects 
in the azimuthal angle correlations and their dependence 
on collision energy. Our conclusions are summarized in 
Section IV. 



II. REMINDER OF THE PHSD APPROACH 

Here we analyze the dynamics of partons, hadrons and 
strings in relativistic nucleus-nucleus collisions within 
the Parton-Hadron-String Dynamics approach [3~l| . In 
this transport approach the partonic dynamics is based 
on Kadanoff-Baym equations for Green functions with 
self-energies from the Dynamical QuasiParticlc Model 



(DQPM) [H, HU which describes QCD properties in 
terms of 'resummed' single-particle Green functions. In 
Ref. [HI, the actual three DQPM parameters for the 
temperature-dependent effective coupling were fitted to 
the recent lattice QCD results of Ref. [35[. The latter 
leads to a critical temperature T c m 160 McV which cor- 
responds to a critical energy density of e c « 0.5 GeV/fm 3 . 
In PHSD the parton spectral functions pj (j = q, q, g) are 
no longer 8— functions in the invariant mass squared as 
in conventional cascade or transport models but depend 
on the parton mass and width parameters which were 
fixed by fitting the lattice QCD results from Ref. (35| . 
We recall that the DQPM allows one to extract a po- 
tential energy density V p from the space-like part of the 
energy-momentum tensor as a function of the scalar par- 
ton density p s . Derivatives of V p w.r.t. p s then define a 
scalar mean-field potential U s (p s ) which enters into the 
equation of motion for the dynamic partonic quasipar- 
ticles. Thus, one should avoid large local fluctuations 
in the potential V p which indeed is solved in the paral- 
lel ensemble method by averaging the mean-field one of 
many events. In the present study we modify the default 
PHSD approach by evaluating the electromagnetic fields 
for each event without averaging the charge currents over 
many (parallel) events. 

The transition from partonic to hadronic d.o.f. 
(and vice versa) is described by covariant transi- 
tion rates for the fusion of quark-antiquark pairs or 
three quarks (antiquarks), respectively, obeying flavor 
current-conservation, color neutrality as well as energy- 
momentum conservation [UHJ- Since close to the phase 
transition the dynamical quarks and antiquarks become 
very massive, the formed resonant 'prchadronic' color- 
dipole states (qq or qqq) are of high invariant mass, too, 
and sequentially decay to the groundstatc meson and 
baryon octets increasing the total entropy. 

On the hadronic side PHSD includes explicitly the 
baryon octet and decouplet, the 0~- and l~-meson 
nonets as well as selected higher resonances as in the 
Hadron-String-Dynamics (HSD) approach J3EH3|- Note 
that PHSD and HSD merge at low energy density, 
in particular below the critical energy density e c rj 
0.5 GeV/fm 3 . 

The PHSD approach has been applied to nucleus- 
nucleus collisions from ^/s^n ~ 5 to 200 GeV in 
Refs. [HI, HH in order to explore the space-time regions 
of 'partonic matter'. It was found that even central colli- 
sions at the top-SPS energy of ^/snn = 17.3 GeV show a 
large fraction of nonpartonic, i.e. hadronic or string- like 
matter, which can be viewed as a hadronic corona. This 
finding implies that neither hadronic nor only partonic 
'models' can be employed to extract physical conclusions 
in comparing model results with data. All these previous 
findings provide promising perspectives to use PHSD in 
the whole range from about tJsnn = 5 to 200 GeV for 
a systematic study of azimuthal asymmetries of hadrons 
produced in relativistic nucleus-nucleus collisions. This 
expectation has been realized, in particular, in the sue- 
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cessful description of various flow harmonics in the tran- 
sient energy range [H, HI] . 




FIG. 1: (Color online) The transverse plane of a non-central 
heavy-ion collision. The impact parameter of the collision 
is denoted by b. The origin O (corresponding to x = 0) is 
referred to as the central point of the maximum overlap. 

The collision geometry for a peripheral collision is dis- 
played in Fig. Q] in the transverse (x — y) plane. The 
reaction plane is defined as the (z — x) plane. The over- 
lapping strongly interacting region (participants) has an 
"almond" -like shape. The nuclear region outside this 'al- 
mond' corresponds to spectator matter which is the dom- 
inant source of the electromagnetic field at the very be- 
ginning of the nuclear collision. Note that in the PHSD 
approach the particles arc subdivided into target and pro- 
jectile spectators and participants not geometrically but 
dynamically: spectators are nucleons which suffered yet 
no hard collision. 

As in Refs. [1, 0> HI, Hll the electric and magnetic fields 
at the relative position R n = r— r„ arc calculated accord- 
ing to the retarded (t n = t — |r — r„|) Licnard-Wiechert 



equations for a charge moving with velocity v 

n v ' 

.BM) - a^Z T J£ (W) (2) 

where the summation runs over all charged quasiparticlcs 
(Z n is the charge number) in the system and a = e 2 /Air = 
1/137 is the electromagnetic constant. Eqs. (U),© have 
singularities for R n = and in the calculations we limit 
ourselves by the condition R n <0.3 fm. 



III. SOURCES OF BACKGROUND 
FLUCTUATIONS 

A. Fluctuations in the proton spectator positions 

Let us consider geometrical fluctuations in the elec- 
tromagnetic field taking into account fluctuations in the 
position of spectator protons. The retarded electric and 
magnetic field evaluated according to Eqs. ([I]),© are pre- 
sented in Fig. [5] for off-central Au+Au collisions at the 
collision energy of ^s^n =200 GeV. The PHSD results 
(including contributions of all quasiparticles) are given 
for the time of the maximal overlap of the compressed 
colliding nuclei which corresponds to t ~0.05 fm/c. As 
noted above the main contribution is coming from specta- 
tor protons. In peripheral collisions the average magnetic 
component orthogonal to the reaction plane < B y > is 
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FIG. 2: (Color online) Probability densities for electromagnetic fields at the central point of off-central Au+Au collisions for 
v / ijviv =200 GeV at the time of maximal overlap as emerging from event-by-event calculations. The results are given for the 
PHSD model (left panel) and for the schematic model with infinitely thin nuclei as in Ref. [2i| (right panel). 
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dominant. The dimensionlcss field magnitude e < B y > 
/m^ 2 ~5 is consistent with earlier calculations within 
the hadronic dynamics of Ultrarelativistic Quantum- 
Molecular Dynamics (UrQMD) [f| and Hadron- String- 
Dynamics (HSD) models as well as with recent cal- 
culations within the partonic HIJING model [3(J. If one 
looks at the field variance (left panel of Fig. [5]) the full 
width of the E y , E x , B x distributions is about a ~ 2/m 2 
for all transverse field components being consistent with 
Ref. [3(|. Here, additional results are plotted also for the 
restricted case when the electromagnetic field is averaged 
over all events in the parallel ensemble as explained in the 
previous section. This procedure has been used before in 
Ref. Q. As seen from Fig. 2 this leads to a suppression 
of the variance for all field distributions by a factor of 
about 3. 

In the right panel of Fig. [5] we mimic results of the 
schematic model in Ref. [29| considering a nuclear col- 
liding system at the time of the maximal overlap as an 
infinitely thin disk. This was simulated numerically by 
an artificial shift of the position of the longitudinal com- 
ponents of all protons at this moment to the plane z = 0. 
As is seen in Fig. 2 (r.h.s.) all field distributions indeed 
increase in width by a factor of about two. A direct com- 
parison of our results to those of Ref. [2^] gives a factor 
of three or even more. This finding completely coincides 
with the results of Ref. as to both the value of the 
width and its origin. 

The estimated strength of the electromagnetic fields 
provides no information about their action on the quasi- 
particlc transport. Let us look at the early time dynamics 
in more detail and introduce a momentum increment Ap 
as a sum of the mean particular increases of the quasi- 
particle momentum dp due to the action of the electric 
and magnetic forces, 



F em = eE + (e/c) v x B, 



(3) 



during the short time interval at the expense of the given 
source, 



Ap(t) = £<dp(t0). 



(4) 



In Fig. [3] the average momentum change of forward mov- 
ing quarks (p z > 0) is shown for three components of 
the electromagnetic force at ^/snn =200 GcV. Note the 
different scales for the solid lines in Fig. [3] that give the 
net momentum change at this energy. It is a remarkable 
fact that the transverse electric and magnetic compo- 
nents compensate each other almost completely. 

Two remarks are in order: First, due to the linearity of 
the electromagnetic force (|3]) with respect to the electric 
and magnetic field, one should not expect a difference 
in quark transport calculations with and without tak- 
ing into account electromagnetic field fluctuations. This 
was demonstrated for quasiparticles earlier in terms of 
the HSD model [27[. Second, if transverse fluctuations 
are characterized by the average strength of the fields, 
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FIG. 3: (Color online) Time dependence of the momentum 
increment of forward moving (p z >0) partons due to the elec- 
tromagnetic field created in Au+Au (y / ijvjv=200 GeV) colli- 
sions with the impact parameter & =10 fm. 



(I-E^yl) and (\B x , y \), certain equalities between compo- 
nents like (\E X \) w (\E y \) ss (\B X \) - as numerically ob- 
tained in Ref. [29| and confirmed in Ref. [3(| - imply 
that similar equalities should hold for the fluctuations. 
Indeed, similar relations follow from our PHSD calcula- 
tions, see Fig. [5] (left panel) where the increment func- 
tions for appropriate field components practically coin- 
cide. We emphasize again that the PHSD transverse 
field components arc not only of comparable strength 
but their action on the quarks approximately compensate 
each other. We recall that in terms of the HSD transport 
model a similar effect has been observed in Ref. [27j ■ For 
a quasiparticle moving along the trajectory x = x{t), this 
compensation in a simplified ID case can be illustrated 
by a short calculation as 



dA 



dAdx 



eE = — e— — ~ — e— — 



Of 



dx dt 



-eBv 



(5) 



i.e. the action of the electric and magnetic transverse 
components is roughly equal and directed oppositely. 

The important advantage of the PHSD approach rela- 
tive to hadron-string models is the inclusion of partonic 
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FIG. 4: (Color online) Time dependence of the momentum increment from 'electric' E c and 'magnetic' B c (two top panels) 
partonic field components, the scalar field and total momentum increment (two bottom panels) for forward moving (p z >0) 
positively charged quarks and p > as a function of time t. The system is Au+Au (at ^/sjvjv=200 GeV) for the impact 
parameter b =10 fm. 



degrees of freedom. In particular, the involved partonic 
fields (of scalar and vector type) showed up to be essen- 
tial to describe the elliptic flow excitation function from 
lower SPS to top RHIC energies [H, and to be a key 
quantity in analyzing the CME. 

The evolution of momentum increments for partonic 
forces is presented in Fig. 2] for off-central Au+Au colli- 
sions at y/SNN =200 GeV. It is seen that the transverse 
'electric' E c and 'magnetic' B c components (two upper 
panels in Fig. @| almost compensate each other. The z 
component is practically vanishing and for t > 10 fm/c all 
quark increments stay roughly constant, i.e. the quark 
phase ends here. The final action of the partonic forces 
is defined by the sum of the forces (right bottom panel) 
which is dominated by the scalar one. 

Apart from the average forces (momentum increments) 
the fluctuations of the forces are of further interest. As 
seen from Fig. [5] the distribution in the quark momen- 
tum deviation Sp = p — (p) in case of scalar forces is well 
collimated with respect to the average trajectory (p) pre- 
sented in Fig. @] but its width increases by about a factor 
of three when proceeding from t =0.05 to 3.0 fm/c (upper 



panel of Fig. [5J- This spread is slightly larger in the x 
component since the derivatives of the scalar mean-field 
are higher in x- than in y-direction. The influence of the 
electromagnetic force on quarks and charged pions is vis- 
ible more clearly (middle and bottom panel of Fig. [5]) 
in the early time corresponding to the maximal overlap 
of the colliding nuclei (t =0.05 fm/c) when the created 
electromagnetic field is maximal. Here, the (8p x ) compo- 
nent is shifted for quarks (middle panels) and even more 
for mesons (bottom panels). This shift decreases in time 
and disappears for t =3 fm/c; at this time the deviation 
distributions for all three components of the electromag- 
netic force are close to a delta function. 

Some general considerations on parity violation in 
heavy-ion reactions are in order here: Since the mag- 
netic field is odd under time reversal (or equivalcntly, 
under the combined charge conjugation and parity CV 
transformation), the time reversal symmetry of a quan- 
tum system is broken in the presence of an external mag- 
netic field. A magnetic field B can also combine with 
an electric field E to form the Lorentz invariant (E • B) 
which changes the sign under a parity transformation. In 
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FIG. 5: (Color online) Probability distribution of the momentum increment <5p = p (p) during the time step St for forward 
moving positively charged quasiparticles at times t = 0.05, 0.1 and 3 fm/c. The distribution emerging from the quark scalar 
potential is shown in the upper panels whereas the distribution stemming from the EM field of quarks and mesons are displayed 
in the middle and bottom panels, respectively. The calculations have been performed for off-central Au+Au collisions at 
y?SNN =200 GeV and impact parameter b =10 fm in the PHSD model. 



the normal QCD vacuum with its spontaneously broken 
chiral symmetry the leading interaction involves the in- 
variant (E ■ B) which enters e.g. into the matrix element 
that mediates the two-photon decay of the neutral pseu- 
doscalar mesons. In the deconfined chirally symmetric 
phase of QCD, the leading interaction term is propor- 
tional to aa s (E-B)(E a -B a ), where E a and B a denote the 
chromoelectric and chromomagnetic fields, respectively, 
and a s is the strong QCD coupling. Both interactions 
are closely related to the electromagnetic axial anomaly, 
which in turn relates the divergence of the isovector axial 
current to the pseudoscalar invariant of the electromag- 
netic field (see Ref. [4(| )■ The evolution of the electro- 
magnetic invariant E • B is shown in Fig. |51 The case of 
Au+Au (y'spf jv=200 GeV) collisions at impact param- 
eter b =10 fm is considered. As seen from Fig. 6 the 
electromagnetic invariant (E • B) is non-zero only in the 
initial time t <0.5 fm/c where the (E • B) distribution 
is quite irregular and its non-zero values correlate well 
with the location of the nuclear overlap region. For later 
times this electromagnetic invariant vanishes in line with 



the electric field space-time distributions Q. Note that 
the quantities plotted in Fig. [6] are dimensionless and the 
scaling factor m% [GeV 4 ] is quite small. 

One should note that in addition to the strong elec- 
tromagnetic fields P, H[ present in non-central collisions, 
very strong color electric E a and color magnetic H a fields 
are produced in the very beginning of these collisions as 
shown in the non-Abelian field theory [4l|. These fields 
can be characterized by a gluon saturation momentum 
Q s and the time ~ 1/Q S - Both fields are parallel to 
each other and directed along the z-axis. This leads to a 
nonzero topological charge Q ~ (E a -B a ) ^0. Since gauge 
fields with Q y^Q generate chirality, they also can induce 
electromagnetic currents along a magnetic field [42| re- 
sulting in the CME. Though a large amount of topologi- 
cal charge might be produced through the mechanism of 
sphalcron transitions, the primary mechanism for topo- 
logical charge Q generation at the early stage is by fluc- 
tuations of color electric and magnetic fields. The decay 
of these fields is essentially governed by the non-Abelian 
dynamics of the Glasma pll |43| which ultimately pro- 
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FIG. 6: (Color online) Space-time evolution of the scalar product of electric and magnetic fields (E ■ B) for Au+Au reactions 
at impact parameter b —10 fm and ^/sjvjv=200 GeV. Note the different scales along the z-axis. 



duces the QGP (close to equilibrium). Unfortunately, 
this possible mechanism for the CME is beyond the po- 
tential of the PHSD model used. We thus may speculate 
about but not prove this mechanism. 



B. Fluctuations in the position of participant 
nucleons 



the collective flows are measured with respect to a third 
plane, the so-called event plane defined by observable 
charged participants in momentum space through the 
harmonic/multipole analysis. More precisely, the flow 
coefficients v n are defined as the n-th Fourier harmonic 
of the particle momentum distribution with respect to 



As noted above (see Fig. [TJ, the interaction region af- 
ter averaging over many events has an almond-like shape; 
the averaged spatial initial asymmetry of the participant 
matter is symmetric with respect to the reaction plane. 
Actual collision profiles, however, are not smooth and 
the symmetry axis in an individual event is tilted due 
to fluctuations (cf. Fig. [7|). The geometry fluctuations 
in the location of the participant nucleons lead to fluc- 
tuations of the participant plane (PP) from one event to 
another, rendering larger coordinate space eccentricities 
which due to pressure gradients are translated into ellip- 
tic flow for the final state particles. Thus, the system of 
the elliptical almond-like shape expands predominantly 
along the minor axis. 

Depending on the location of the participant nucle- 
ons in the colliding nuclei at the time of the collision, 
the actual shape of the overlap area may vary. As is 
seen from Fig. [7J due to fluctuations the overlap area 
in a single event can have, for example, a rotated tri- 
angular rather than an almond shape. Note that an al- 
mond shape is regained by averaging over many events 
for the same impact parameter. However, in experiment 
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FIG. 7: (Color online) Projection of a single peripheral 
Au+Au (200 GeV) collision on the transverse plane. Specta- 
tor and participant nucleons are plotted by empty and filled 
circles, respectively. The reaction plane (RP) projection cor- 
responds to x axis. Transverse axes of the participant plane 
(PP) are marked by stars (x*,y*). 
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the particular momentum event plane "F^, 

(v n ) = (cos[n(4> - *„)]) (6) 

where tp — arctan(p y / p x ) is the azimuthal angle of the 
particle momentum p in the CM frame and angular 
brackets denote a statistical average over many events. 




-0.8 -0.6 -0.4 -0.2 -0.0 0.2 0.4 0.6 0.8 

nw 1% 

1 n 

FIG. 8: (Color online) Distribution in the event plane angle 
for different harmonics ty n calculated with retarded magnetic 
and electric fields. Grey histograms show the results for re- 
spective calculations without fields. 

One should note that all azimuthal correlations are 
not only due to the collective flow. The early-time two- 
particle spatial correlations probe both the event geom- 
etry (fluctuating in individual events) and genuine lo- 
cal pair correlations referred to as 'nonflow' correlations. 
The Fourier decomposition Q is not enough to disentan- 
gle these two- contributions. A possible solution of the 
connection between flow fluctuations and initial state cor- 
relations is given by the cumulant expansion method |44j 
using two- and four-particle correlation measurements of 
the harmonic flow coefficients. However, this method is 
beyond the scope of the present study. 

The distributions in the event plane angle for different 
harmonics are shown in Fig. [5] for the freeze-out case. 
All distributions are symmetric with respect to the point 
<F„ = which corresponds to the true reaction plane. As 
is seen, the event plane angle ty n determined from the 
"n-th" harmonic is in the range < < 27r/n and 



fluctuations of several lowest order harmonics have com- 
parable magnitudes. Inside this region >Fi has two max- 
ima at ^ n =0 and 7r corresponding to forward-backward 
emission. The even components "I^, ^4 have a rather 
prominent maximum for ^„ = indicating the local na- 
ture of fluctuations, but the odd harmonics ^I^,^ are 
practically flat. This may be easily understood since the 
odd moments of the spatial anisotropy purely originate 
from fluctuations while the even ones are combined ef- 
fects of fluctuations and geometry. As a consequence, 
if one defines the spatial anisotropy parameters with re- 
spect to the pre-determined reaction plane, the event- 
averaged ones vanish for all odd moments but not for the 
even moments. 

The histograms in Fig. [8] are calculated from a sample 
of 3.10 events taking into account magnetic and electric 
field fluctuations. Similar calculations without fields arc 
shown in the same figure by the gray histograms which 
are hard to distinguish from the previous ones. In other 
words, there is no additional "tilting" effect by electro- 
magnetic fields as expected in Ref. [2^, [3(| ■ This is due 
to the compensation of the transverse electromagnetic 
components as explained above. 



C. Two-particle angular correlations 

An experimental signal of the local spontaneous parity 
violation is a charged particle separation with respect to 
the reaction plane j45|. It is characterized by the two- 
body correlator in the azimuthal angles, 

jij = (cos(V>i + tpj - 2^ RP )) (7) 

= (cOS(V>j - *B,p) COS(V>j - $rp)) 

- (sin(V> 4 - ^ RP ) sin(-0j - *fip)) 

where ^ pp is the azimuthal angle of the reaction plane 
defined by the beam axis and the line joining the centers 
of the colliding nuclei and subscripts of 7y represent the 
signs of electric charges being positive or negative. The 
averaging in Eq. ([7]) is carried out over the whole event 
ensemble and cos and sin terms in ([7]) correspond to out- 
of-planc and in-plane projections of 7^ . 

As was pr oposed inRefs. [13, El and more elaborated 
in Ref. [22j , a possible source of azimuthal correlations 
among participants is the conservation of the transverse 
momentum which might give rise to a contribution com- 
parable with the measured CME. Transverse momentum 
conservation (TMC) introduces back-to-back correlations 
for particle pairs because they tend to balance each other 
in transverse momentum space. A large multiplicity of 
particles will dilute the effect of these two-particle corre- 
lation. Furthermore, this correlation should be stronger 
in plane than out of plane due to the presence of the ellip- 
tic flow. Nevertheless, the TMC provides a background 
for the CME that should be properly quantified. 

From quite general considerations - making use of the 
central limit theorem and describing particles thermody- 
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namically - one can derive the following simple expression 
for the two-particle correlator [22j 



jij = (cos(ipi + Ipj)) 



V 2 (j> 



t/ acc 



N (p? 



>full 



(8) 



where N = N + + Nq + N- is the total number of all pro- 
duced particles (in full phase space). At y / s]vjv =200 
GeV it can be approximated as N w (3/2) N c h ~ 
21 N part where N part is calculated dynamically in 
our model as well as the momentum-dependent factors 
for full phase space and the ratio to the measured ac- 
cepted phase space. It is of interest to note that the 
proportionality of the CME to the elliptic flow v-i seen in 
Eq. follows also from more elaborated consideration. 
In particular, the chiral magnetic effect in the hydrody- 
namic approach and in terms of a holographic gravity 
dual model (see Ref. [46|) predicts a linear dependence 
of the CME on V2 with more sophisticated coefficients 
which depend on the axial anomaly coefficient and the 
axial chemical potential as well as on dynamics of fluids 
through the particle density, baryon chemical potential 
and pressure. 

Experimentally, the same-sign correlator 7 SS is defined 

as the average of 7++ and 7 by assuming that the 

momentum balance is shared equally among the charges 



1 



7s 



(7++ +7—) 



v 2 (Pt)l cc 

'N (P 2 t )full 



(9) 



In practice, only a subset of particles is measured. In this 
case some of the momentum balance stems from unmea- 
sured particles and one might expect 7 SS <C V2/N [lj|. 
In the STAR experiment [47| tracks were measured for 
the central two units of rapidity. However, the initial 
colliding beams approached with ±5.5 units of rapidity 
and more than 50% of the charged particles tracks have 
rapidities outside the STAR acceptance. These particles 
can serve as a source of momentum, which may quench 
the momentum conservation condition thus reducing the 
magnitude of j ss . However, the transverse momentum of 
a given track is more likely to be balanced by neighboring 
particles, which have similar rapidities. This is particu- 
larly true when considering the components of the mo- 
menta responsible for elliptic flow. We conclude that this 
effect should be more essential for lower collision energy. 

The direct comparison of the momentum conservation 
effect ([H]) on the CME observable is presented in Fig. [5] 
for the top RHIC energy. The total (rather than trans- 
verse) momentum conservation is inherent in the PHSD 
model. In the actual calculations the experimental ac- 
ceptance p t >200 MeV/c is taken into account; as seen 
from the upper part of Fig. [9] the ccntrality dependence 
of the elliptic flow (^2) for charged particles is rather well 
reproduced by PHSD. However, the experimental same- 
sign correlator 7 SS is underestimated substantially. We 
note that the experimental acceptance essentially influ- 
ences the momentum-dependent ratio (pt)acc/ (Pt) full- ^ n 
reality the difference in 7 SS should be even larger as dis- 
cussed above. This point is in agreement with the full 
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FIG. 9: (Color online) PHSD centrality dependence of the 
elliptic flow (top part) and angular correlators 7 SS and S S3 
of charged particles (lower part) from Au+Au at V / Sjvjv=200 
(middle) and 7 (bottom) GeV from the transverse momentum 
conservation according to Eqs. ^ and (|10[) . The experimen- 
tal data points for V2 and 7 SS , 8 SS are from Refs. [48|] and 
[13], respectively. 



HSD calculation of the hadronic background within the 
CME studies in Ref. [27j]. 

A similar analysis for the lower energy ^snn =7.7 
GeV is presented in the lower panel of Fig. [5] Unfortu- 
nately, measured data for the centrality dependence of 
the elliptic flow are not available at this energy but the 
PHSD calculated average (1*2) for minimum bias colli- 
sions is only slightly below the experiment [39| due to ne- 
glecting a baryon mean-field potential (see also the end of 
subsection IIIID[) . The calculated correlation 7^ strongly 
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differs from the measured values having even the oppo- 
site sign! One should note that in this case the same- and 
opposite sign components are almost equal to each other 
(i.e. there is no charge separation effect). This observa- 
tion is also nicely reproduced within the HSD model at 
this energy (c/. Ref. 

It is of further interest to consider the average cosine 
of the transverse angle difference which is independent of 
the reaction plane 



% = (cos(ipi - ijjj)) = - 



(pi 



N(p 2 t 



full 



(10) 



where the last equality is obtained from the transverse 
momentum conservation p2| . As follows from the com- 
parison between (J5J and (|TU)) . the correlator <5„ differs 
from only by the elliptic flow coefficient V2 and is ex- 
pected to be more sensitive to the TMC. As one can 
see from Fig. [5] this estimate of S ss is too large and 
hardly consistent with appropriate experimental data 
from Fig. HU 

Thus, the considered angular correlation 7„ s is gener- 
ated by a combination of momentum conservation, which 
causes particles to be preferably generated in the oppo- 
site direction, and elliptic flow which gives more particles 
in the ±x direction than in the ±y direction. However, 
this source is by far not able to explain the observed 
pion asymmetry in the angular correlation. In addition, 
the considered TMC is blind to the particle charge and 
cannot disentangle same-sign and opposite-sign pair cor- 
relations. 



D. Electric charge fluctuations in the transient 
stage 

The almond-like fireball created in the early collision 
phase then expands in an anisotropic way, however, the 
spatial anisotropy is reduced with increasing time. In 
this transient stage the electromagnetic field is strongly 
reduced since the spectator matter is flying away from 
the formed fireball. The pressure gradients act predomi- 
nantly in the reaction plane resulting in elliptic flow V2- 

Strong interactions in this phase might produce signifi- 
cant fluctuations in energy density (temperature) , trans- 
verse momentum, multiplicity and conserved quantities 
such as the net charge. In the plasma phase in a mag- 
netic field an electric quadrupole can be formed due 
to chiral anomaly and as a signal of that the elliptic 
flow difference between 7r + and 7r~ mesons is predicted 
V2(tt + ) < V2(n~) Certainly, to check that the in- 

fluence of hadronic transport on observablcs should be 
taken into account. 

Furthermore, the CME [l[ predicts that in the presence 
of a strong electromagnetic field at the early stage of 
the collision the sphaleron transitions in a hot and dense 
QCD matter induce a separation of charges along the 
direction of the magnetic field which is perpendicular to 
the reaction (x — z) plane. This charge separation results 



in the formation of an electric dipole in momentum space 
which breaks parity. Being interested essentially in the 
quark phase, we investigate in this subsection to what 
extent such an electric dipole can be generated by back- 
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FIG. 10: (Color online) Probability distribution in the mag- 
nitude Qd of the generated electric dipole at freeze-out for all 
partons. The top and bottom panels correspond to calcula- 
tions without and with the electromagnetic field, respectively. 
The system is Pb + Pb at ^/sjvjv = 200 GeV at impact pa- 
rameter b— 10 fm. 

ground statistical and electromagnetic field fluctuations. 

Let us quantify the dipole defining the plane Q c i of the 
quark distribution in the transverse momentum space. 
The magnitude Q c \ and azimuthal angle VPd of this vec- 
tor can be determined in a given event as follows: 



Qd cos * c i = ^ qi cos ipi , 

i 

Qd sin * c i = 2J Qi sin ipi 



(11) 



where the summation runs over all charged particles in 
the event with the electric charge qt and azimuthal an- 
gle t/ji of each particle. Note that Eq. (fTTj) describes 
the dipole shape of charged particles (quarks or hadrons) 
without any reference to the charge separation. 
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As seen from Fig. [KjJ the average charge magnitude 
of the electric dipole Q c \ at the moment of the maxi- 
mum nuclear overlap (t =0.05 fm/c) is about 4 charge 
units with dispersion ctqx ~2. At this moment the sys- 
tem is in the quark phase having on average an almond- 
like shape and therefore is expanding preferentially along 
the x— axis. Note that according to Eq. (fTTj) quark net 
electric charges with \qi\ <1 are considered. Thus, the 
number of quarks involved in the dipole is large. In the 
next step (t = 0.1 fm/c) of the expansion stage (dot- 
dashed line in Fig. [TU]) the Q c i distribution is getting 



AuAu, \fS^=200 GeV, b=10 fm 




FIG. 11: (Color online) Electric dipole evolution of charged 
partons in Q c i — t presentation for Au+Au collisions at 
^Jsnn = 200 GeV with taking into account the electromag- 
netic field. 

broader with a noticeable increase of Q c i- At t =10 fm/c 
the quark-gluon phase transforms predominantly into the 
hadronic phase through the dynamical coalescence mech- 
anism and the Q c \ distribution becomes narrower again. 
The influence of the electromagnetic field on this evolu- 
tion is very weak (compare the top and bottom panels in 

Fig. nnj). 

The whole evolution of the electric dipole is seen more 
clearly in the 3D representation in Fig. QTJ Indeed, the 
Qd distribution has a pronounced peak formed shortly 
after the collision, then the magnitude of the Q c \ charge 
distribution is minimal during about 3 fm/c to testify 
that a large-in-charge subsystem is formed. After that 
the Q c i-distribution is getting narrower because during 
the expansion the Qci value slightly decreases due to par- 
ton hadronization; the maximum of the probability dis- 
tribution increases and then stabilizes after t salO fm/c. 

Similarly to the dipole, one can define a charged 
quadrupole formed in heavy-ion collisions as follows: 

Q c2 cos 2* c2 = 2J It cos , 

i 

Q c2 sin 2* c2 = y^gj sin 2ipi . (12) 

i 

Characteristics of the time evolution of the electric 
dipole and quadrupole (formed in semi-central Au+Au at 
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FIG. 12: (Color online) Distribution in the magnitude of the 
charge dipole (upper panel) and quadrupole (bottom panel) 
calculated for charged particles at times t= 0.05, 1, 10 and 
40 fm/c. The system is Au + Au at ^snn— 200 GeV and 
impact parameter b= 10 fm. 



y?SNN = 200 GeV collisions) for all (parton and hadrons) 
charged quasiparticles in the midrapidity range are pre- 
sented in the next two figures. The Q c \ minimum ob- 
served in the pure partonic case (see Figs. [TU] and fTTj) 
survives also in this case. The magnitude of the n =2 
quadrupole harmonic presented in Fig.ll2p is close to that 
for the dipole n =1, i.e. Q C 2 ~ Q c i and their maximal 
values extend to values of 30-40. As is seen from Fig. [T3l 
the distributions in the reaction plane angle for the elec- 
tric quadrupole 1 F C 2 are rather flat during the whole evo- 
lution while the electric dipole angle distribution is 
flat only in the partonic phase (see t =0.05 fm/c in Fig. 
1 1 3|) but in the hadronic phase the distribution resembles 
that for the directed flow (c.f. Fig. [5]). The main axis 
of and can randomly be parallel or anti-parallel 
to the minor axis of the almond. Like in the dipole case 
(cf. Fig. [TU)) the electromagnetic field has no sizable in- 
fluence on the characteristics of the electric quadrupole. 
When the collision energy yijvjv decreases the behavior 
of the dipole and quadrupole distributions in the mag- 
nitudes Qd and Q C 2 and the angle practically do not 
change besides some structure in the \l/ c i distribution. 
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FIG. 13: (Color online) The same as in Fig. [12] but for the 
distribution in the event plane angle of the charge dipole (up- 
per panel) and quadrupole (bottom panel) at times t= 0.05, 
1, 10 and 40 fm/c. 



As seen from Fig. [TJ]at ^/JWn =11.5 GeV back-to-back 
correlations - as specific for the directed flow - are man- 
ifested. This is mainly due to the proton contribution 
which becomes noticeable at low collision energy. 

Thus, the statistical fluctuations of "normal" matter in 
the presence of the retarded electromagnetic field do not 
result in a sizeable formation of a deformed subsystem of 
dipole- or quadrupole-shape during the evolution of the 
heavy-ion collision. 

We point out, however, that such subsystems might 
be formed in non-trivial topological systems due to the 
chiral anomaly effect. In particular, it happens when 
a quark experiences both a strong magnetic field and a 
top olog ically nontrivial gluonic field such as an instan- 
ton [50J ■ The inherent asymmetry - when both instanton 
and magnetic field are present - can lead to the devel- 
opment of an electric dipole moment. Physically, it can 
be understood as the result of two competing effects: the 
spin projection produced by a magnetic field and the chi- 
rality projection produced by an instanton field. Such a 
consideration is beyond the scope of our present micro- 
scopic study. 

It is of interest that the axial anomaly in a strong 
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FIG. 14: (Color online) Angular distribution of charged parti- 
cles for the charge quadrupole (solid line) and dipole (dashed) 
subsystems calculated for Au+Au (^/snn = 11.5 GeV) colli- 
sions. 



external magnetic field induces not only the CME but 
also the separation of the chiral charge. The coupling 
of the density waves of electric and chiral charge results 
in the Chiral Magnetic Wave and can induce a static 
quadrupole moment of the electric charge density p9| . 
This chiral magnetic wave results in the degeneracy be- 
tween the elliptic flows of positive and negative pions 
leading to t>2(7r~) > V2(tt + ), which was estimated on 
the level of ~ 30% for mid-central Au+Au collisions at 
y/sNN =11 GeV [49|. Our PHSD calculations give about 
6% which is quite comparable with the recently mea- 
sured value of 10% [5l[ and essentially smaller the pre- 
diction of Ref. . Noteworthy that the v 2 degeneracy 
in the PHSD version used is only due to different clastic 
and inelastic cross sections for tt + and 7r~ mesons but 
without taking into consideration the (small) mean-field 
pion-nucleus potential. The elliptic flow analysis of the 
difference between particles and anti-particles (including 
kaons and baryons alongside with pions) shows that this 
difference is coming mainly from the hadronic mean-field 
potential 



52|. 



Recently these v-i data have been also 
successfully explained in terms of a hybrid model, which 
combines the fluid dynamics of a fireball evolution with 
a transport treatment of the initial and final hadronic 
states [531 ] . Therefore, there is not much room for the 
contribution from a transient charged quadrupole due to 
the chiral magnetic wave. 



E. Charge balance functions 

In the formation of the charged dipole and quadrupole 
there is no information about a possible charge separation 
which could result in an electric driving force. In prin- 
ciple such information can be provided by the balance 
function which is based on the idea that charge is locally 
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conserved when particles are produced pair-wise. In the 
subsequent expansion of the system and rescattering of 
the charge carriers, which in principle can be hadronic 
or partonic, the balancing partners are then spread out 
within some finite distance to each other. The original 
correlation in space-time transforms into a correlation in 
momentum space in the final hadronic emission profile. 
Therefore, the motion of the balancing partners suffers 
from the collective expansion of the system and diffusion 
due to the collisions with other particles. The study of 
charge-balance correlations hence gives insight into the 
production and diffusion of charge. In particular, it is 
expected that the balance function is sensitive to the de- 
layed transition of the quark-gluon phase to a hadronic 
phase [EH- 

Thus, whenever a positive charge is created, a negative 
charge arises from the same point in space-time and both 
particles then tend to be focused in the same rapidity and 
azimuthal angle by collective flow. This results in a cor- 
relation between positive and negative charges, i.e., for 
every positively charged particle emitted at an angle 
there tends to be a negatively charged particle emitted 
with ip~ w and similar rapidity. Charge balance func- 
tions [HI represent a measure of such correlations, and 
have already been investigated as a function of relative 
rapidity for identified particles and for relative pseudo- 
rapidity ?y for non-identified particles [13, [H[ . Generally, 
the balance function B(p a \pb) is a six-dimensional func- 
tion of the particle momenta. In the context of studies 
of the separation of balancing charges, the discussion is 
reduced to the difference (pi — P2). In particular we will 
focus on the charge balance function in relative pseudo- 
rapidity <5?y, i.e. B(p a \pb) — > B(5r],r] w ) and similarly for 
the azimuthal angle i[> [561 ] . 

Charge balance functions are constructed in such a way 
that like-sign subtractions statistically isolate the charge 
balancing partners, 



B(Sr],r) w ) 



If N + -(8r),r) w ) - N ++ (5rj,T] w ) 



2 V N + 
N- + (6r),r] w ) - N — {Srj,T] w ) 
A_ 



, (13) 



where the conditional probability (5ri,r) w ) counts 

pairs with opposite charge which satisfy the criteria 
that their relative pseudorapidity Sr/ = r/ + — ?y_ in a 
given pseudorapidity window is rj w , (Srj <G 77^), whereas 
jV + (iV_) is the number of positive (negative) particles 
in the same interval. Here the angular bracket repre- 
sents averaging over the events. Similarly for N ++ ,N 

and N |_. The factor 1/2 ensures the normalization of 

B(Sr],r] w ). All terms in Eq. (|13|) arc calculated within 
PHSD using pairs from a given event and the resulting 
distributions are summed over all events. 

Both balance-function and charge-fluctuation observ- 
ables are generated from one-body and two-body observ- 
ables which necessitates that they may be expressed in 
terms of spectra and two-particle correlation functions. 



The charge fluctuation is a global measure of the charge 
correlation and the balance function is a differential mea- 
sure of the charge correlation; it therefore carries more 
information. Writing N± = (N±) nw +SN±, where (.■■) r)w 
denotes the average in the phase-space region rj w , it is 
easy to show that 



((Qch - (Qch)f 



eh I 



dSr) B(Sr)\r) u 



(14) 



where Q ch = N + - N- and N ch = A + + A_ . 

From this example, one can readily understand how 
balance functions identify balancing charges. For any 
positive charge, there exists only one negatively charged 
particle whose negative charge originates from the point 
at which the positive charge was created. By subtract- 
ing from the numerator the same object created with 
positive-positive pairs, one is effectively subtracting the 
uncorrelated contribution from the distribution. 

It is expected that charge balance functions are sensi- 
tive to the separation of balancing charges in momentum 
space and give insight into the dynamics of hadroniza- 
tion [54j]. Indeed, such a pair is composed of a positive 
and negative particle (or particle and antiparticlc) whose 
charge originates from the same point in space-time. If a 
quark-gluon plasma results in a large production of new 
charges late in the reaction, a tight correlation between 
the balancing charge-antichargc pairs would provide evi- 
dence for the creation of a novel state of matter. 

The time evolution of the Srj- and <50-dcpendent charge 
balance function is demonstrated in FigHS] for central 
(left panels) and peripheral (right panels) Au+Au col- 
lisions at y/sNN = 200 GeV. The times t =0.05 and 
1 fm/c correspond to the developed quark phase which 
ends at about t =10 fm/c (cf. Fig. i| while at t =40 
fm/c the system is in a purely hadronic phase. We re- 
call that hadronization in the PHSD model is realized 
via a crossover transition and quasiparticle rescattering 
is included in both the partonic and hadronic phase. 

As follows from Fig. [15] there is a very small difference 
in the time evolution of charge balance function B(Sjj) 
at both ccntralitics. As to B(S(j>) a rather clear enhance- 
ment is seen for central collisions at Sr] ~ 8(f> ~0 while 
this dependence essentially weaker for peripheral colli- 
sions. This observation is in qualitative agreement with 
experiment but enhancement effect is too small. Thus 
the expectation of a high sensitivity of the balance func- 
tion to the hadron phase transition and particle diffusion 
seems somewhat too optimistic. 

The direct comparison with experiment of the charge 
balance function is presented in Fig. [16] for Au+Au col- 
lisions at y/sNN =200 GcV calculated within PHSD for 
central and peripheral (b =10 fm) collisions. Here the 
charge conservation law is locally fulfilled in each quasi- 
particle collision. Hadronic resonance decays are taken 
into account in PHSD but corrections due to final state 
interactions for small Stj) are computationally very in- 
volved (and uncertain). 
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8(j) 5<t> 

FIG. 15: (Color online) The balance function for charged quasiparticles (quark-antiquarks and ± pions) with |?)+/-! < 1 from 
central (left panel) and peripheral (right panel) Au+Au (^/snn = 200 GeV) collisions at t =0.05, 1, 10 and 40 fm/c. 



As is seen from Fig. [TO] there is a maximum for 
6r) ~ S<j) ~ but no large difference is observed for dif- 
ferent centralities (apart from the 5<p distribution). This 
finding is in agreement with the kinetic results of UrQMD 
and HIJING model calculations [HH of the width of the 
balance function in terms of 8rj which also show no nar- 
rowing of the peak for central collisions as observed in 
experiments. The inclusion of the electromagnetic field 
(the dotted line in the bottom panel in Fig. [16] for b =10 
fm) practically shows no influence and does not result in 
any additional focusing in central collisions. 

We find that PHSD does not describe quantitatively 
the experimental balance functions. There are some 
claims that the blast-wave model can resolve this dis- 
crepancy [13, EM H3- The blast- wave model is in fact 
a parametrization of the kinetic freeze-out configuration 
motivated by a hydrodynamical model for the system de- 
scribed in local thermal equilibrium. The system is then 
completely characterized by the collective velocity profile, 
freeze-out temperature and the freeze-out surface which 
is usually associated with some volume. Generally, the 
blast-wave model parameters may be varied in a large 
parameter space to fit experimental data. These single- 
particle freeze-o ut p roperties can e.g. be parameterized 
as suggested in |59j to study, for example, the evolution 



of flow. However, the change in the kinetic freeze-out 
temperature and the increase of collective flow alone fail 
to explain the observed focusing of the balance function 
for more central collisions (l8j | . 

With regard to charge-balance correlations, the blast- 
wave model needs additionally to incorporate local charge 
conservation. This can be achieved in the following way: 
Instead of generating a single particle at a time, an en- 
semble of particles with exactly conserved charges is gen- 
erated in such a way to remain unchanged the single- 
particle distributions. For the relative distribution of the 
pairs within an ensemble, a Gaussian distribution is as- 
sumed with dispersions and <t| for rapidity and trans- 
verse angle, respectively. Treating these dispersions as 
free parameters at every centrality it is possible to tune 
the narrowing effect for central events |l8[. It is of in- 



terest that at the exactly central collision = = 
and they strongly grow with impact parameter reach- 
ing cr? «0.6 and cr^/ir ~ 0.4 for centralities about 
70% [18( 1 . However, the additional assumption that bal- 
ancing charges at the freeze-out are strongly correlated in 
momentum space is in conflict with the basic model as- 
sumption on thermodynamic equilibrium of the system. 

In Ref. [l|| the charge balance function was applied for 
the analysis of the CME. The charge separation between 
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FIG. 16: (Color online) The balance function for pseudorapid- 
ity (top) and azimuthal angle (bottom) of charged pions with 
j?7 + /_| < 1 and p t >0.2 GeV/c from central and peripheral 
Au+Au collisions at ^/snn =200 GeV. The dotted lines for 
8(j) correspond to calculations including the electromagnetic 
field effects. The experimental data points are from Ref. fisl |. 



opposite-charge and same-charge two-pion correlators 



was defined as 



7-P 



7P 



2( 2 7+- -74 



7 — ) = 7+- - Iss, (15) 



where angle brackets in ([Jj include the balance function 
B(p + \p_) as a weight factor for the balancing charges. 
The quantity 7p can be estimated from available exper- 
imental data Q. Since the PHSD is not successful in re- 
producing the charge balance function, there is not much 
sense to apply it for the charge separation jp . As follows 
from the comparison between the STAR data and the 
blast-wave model (including correlations and rescaling) 
results in reproducing the experimental normalization; 
the charge balance correlations for the relativistic charge 
separation are of the same size as the experimental sig- 
nal and exhibit a similar qualitative behavior with re- 
spect to the centrality dependence [l||]. The authors of 
Ref. [HI claim that their results are solid on the level of 
10-20%. The calculation of uncertainties originates pre- 
dominantly from the particular paramctrization of both 



the blast-wave model itself and, in particular, of the cen- 
trality dependence of the charge separation B(p a \pb) in 
the azimuthal angle. However, the considered breakup 
physics differs significantly from more realistic scenarios 
as has been shown recently in Ref. [60j; the freeze-out 
temperature and baryon chemical potential - defining the 
chemical composition of the system - noticeably depend 
on centrality. Furthermore, there is some inconsistency 
in using a reaction-plane independent fit of the balance 
function for the CME signal where azimuthal angles are 
measured with respect to the reaction plane. It is also un- 
clear how the extracted parameters change with the col- 
lision energy, however, the first preliminary STAR date 
on the collision-energy dependence of the balance func- 
tion have been published recently 



611. Thus, a further 



careful study of this issue is needed. 



IV. THE CME OBSERVABLE 

The experimental signal of the possible CME is the az- 
imuthal angle correlator calculated according to Eq. ([7J . 
The experimental acceptance |^| < 1 and 0.20 < p t < 2 
GeV has been also incorporated in the theoretical PHSD 
calculations. Note that the theoretical reaction plane is 
fixed exactly by the initial conditions rather than by a 
correlation with a third charged particle as in the exper- 
iment [HI- Thus, within PHSD we calculate the observ- 
able ([7]) as a function of the impact parameter b or the 
centrality of the nuclear collisions which should be con- 
sidered as a background of the CME signal. Comparison 
of measured angular correlator with calculation results 
are presented in Fig. [T7] Calculation of this characteris- 
tic is a very CP time consuming process and the proper 
calculation error bars are plotted in Fig. [TTJ 

At the lowest measured energy ^/s^n =7.7 GeV the re- 
sults for oppositely and same-charged pions are very close 
to each other and show some enhancement in very pe- 
ripheral collisions. The centrality distributions of 7^ are 
reasonably reproduced by the PHSD and HSD calcula- 
tions presented in the same picture. Note that the scalar 
quark potential is not zero at this low energy but absent 
in the HSD model. The striking result is that the case 
of ^snn =7.7 GeV drastically differs from ^/snn =200 
GeV (cf. the right bottom panel in Fig. [17)) . The picture 
quantitatively changes only slightly when one proceeds 
to -y/sjviv =11.5 GeV (see the right top panel in Fig. [T7|) 
though the value of 7^ at the maximum (centrality 70%) 
decreases a little bit in the calculations. Experimental 
points at larger centrality are not available but are of 
great interest. In addition, one may indicate a weak 
charge separation effect in the experimental data because 
statistical error bars are very small (less than the sym- 
bol size). Unfortunately, the calculated error bars are 
rather large to specify the charge separation effect. The 
influence of the electromagnetic field here is negligible. 
The calculated and measured correlation functions for 
oppositely and same charged pions are shown in Fig. [TTJ 
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centrality [ % ] centrality [ % ] 

FIG. 17: (Color online) Angular correlations of opposite- and same-charge pions in azimuthal angle for Au+Au collisions at 
^/s~NN =7.7, 11.5, 39 and 200 GeV as a function of centrality. The full symbols are preliminary STAR data [13] as well as 
published STAR data for ^s]7n = 200 GeV 0. 



for the available three BES energies. The case for the 
top RHIC energy ^snn =200 GeV is also presented for 
comparison. 

If one looks now at the results for -^/sjvjv =39 GeV, 
the measured same- and oppositely charged pion lines 
are clearly separated, being positive for the same-charged 
and negative for the oppositely charged pions to be 
strongly suppressed. The PHSD model is not able to 
describe this picture and overestimates the data with in- 
creasing energy. These growing large values of jij are due 
to the scalar parton potential which increases with the 
collision energy. The HSD version predicts a very small 
effect in qualitative agreement with our earlier analy- 
sis [13 • Though both models provide the charge sep- 
aration essentially smaller than the measured one, the 
PHSD has a satisfying feature: the same-charge points 
are above the oppositely charged ones in agreement with 
experiment. The same situation is observed in the case of 
y/sjVTV =200 GeV; a small difference between them is seen 
in very peripheral collisions: the oppositely charged cor- 
relation jumps to zero at centrality ~ 70% for ^/snn =39 
GeV while corresponding data at 200 GeV are not avail- 
able. 

Though the results at -^/stvjv =7.7 and 11.5 GeV 
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FIG. 18: (Color online) Angular correlations of opposite- and 
same-charge pions for the cosine of the difference in the az- 
imuthal angles for Au+Au collisions at ^snn =200 GeV as 
a function of centrality. The experimental data points are 
from gj]. 



roughly can be considered as a background of the CME, 
at higher energies it is impossible to identify the true ef- 
fect of the local parity violation as the difference between 
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measured and PHSD results. The PHSD model [H[ in- 
cludes directly the dynamics of quark-gluon degrees-of- 
freedom which are becoming more important with in- 
creasing energy. We recall that the growing importance 
of the repulsive partonic mean field - illustrated earlier by 
the rise of the elliptic flow explained convincingly in the 
PHSD model [38, 39[ - results here in an overestimation 
of the CME background. 
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FIG. 19: (Color online) Angular correlations of cos-(out-of- 
plane) and sin(in-plane)-projections for Au+Au collisions at 
y/SNN =200 GeV as a function of centrality. The experimen- 
tal data points are from |J. 

In Fig.[T8]thc results for the average cosine of the differ- 
ence in the azimuthal angles Sij are presented. The mea- 
sured centrality dependence for the same charge pions 
is flat and practically consistent with zero while that for 
unlikely charged particles is a monotonic increasing func- 
tion with impact parameter. However, the PHSD calcu- 
lations clearly overestimate the experimental points 
We note in passing that the PHSD results for Au+Au 
collisions at the energy ,Jsnn= 200 GeV turn out to be 
astonishingly close to the appropriate experimental data 

at y5]vlv =2.76 TeV [HEil. This fact indicates that the 
strength of the repulsive scalar quark potential in PHSD 
might be presently overestimated. 

In accordance with Eq. ([7]), one can separate the in- 
plane and out-of-plane components using experimental 



results for 7^ and 5{j. Such separation together with 
PHSD calculation results is presented in Fig. Q1J] for the 
same charge and opposite charge pions. As was first 
noted in Ref. [l6| and is seen in the upper panel of Fig.lT^l 
for the same-charge pairs the sinus term is essentially zero 
whereas the cosine term is finite. This tells us that the 
observed correlations are actually in-plane rather than 
out-of-plane. This is contrary to the expectation from 
the chiral magnetic effect, which results in same-charge 
correlations out of plane. In addition, since the cosine 
term is negative, the in-plane correlations are stronger for 
back-to-back pairs than for small angle pairs. The PHSD 
does not reproduce these features. We see also that for 
opposite-charge pairs the in-plane and out-of-plane cor- 
relations are virtually identical. As was stated in 
this is difficult to comprehend since there is a sizable el- 
liptic flow in these collisions. Nevertheless, the PHSD 
model predicts very close in-plane and out-of-plane dis- 
tributions for opposite-charge pairs due to scalar parton 
potential and at the same time nicel y rep roduces the var- 
ious harmonics of charged particles [3|| [39| . This feature 
is not reproduced in the HSD. 

We close this Section with some more general remarks. 
As follows from the results presented in Figs. 1171181191 
an additional sizable source of asymmetry is needed for 
both in-plane and out-of-plane components rather than 
only an out-of-plane component as expected from the 
CME. As discussed in the Introduction, the vacuum non- 
trivial topological structure (as a genuine source of the 
CME) leads to the picture of a topological 9- vacuum of 
non-Abelian gauge theories. The 0-term in the QCD La- 
grangian explicitly breaks P and CP symmetries of QCD. 
However, stringent limits on the value of 9 < 3 x 10~ 10 
deduced from the experimental bounds on the electric 
dipole moment of the neutron [(32[ practically indicate the 
absence of global P and CP violation in QCD. Reference 
to the local P- and CP-odd effects due to the topologi- 
cal fluctuations characterized by an effective 9 = 0(x, t) 
varying in space and time f63j does not provide much 
hope. In addition, partons near the phase transition are 
not chiral (as typically assumed) but massive degrees- 
of-freedom in the PHSD in agreement with lattice QCD 
calculations. The finite mass of the partons washes out 
the chirality effect. 



V. SUMMARY AND OUTLOOK 

In this study we have investigated several effects that 
might contribute to the observed chiral magnetic effect 
(CME) in relativistic nucleus-nucleus collisions on the ba- 
sis of event- by- event calculations within the PHSD trans- 
port approach. The individual results can be summarized 
as follows: 

• Our study shows that fluctuations in the position 
of quasiparticles can manifest themselves in differ- 
ent interaction stages and in different ways. Since 
the electromagnetic field generated by spectators is 
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dominant at the early stage, the fluctuation in their 
position results in a noticeable fluctuation in the 
strength of the electromagnetic field. However, the 
fluctuation spread is not so large as expected in 
the estimate from Ref. [1^] and its influence on ob- 
servables is negligible; in particular, the event plane 
angle is not tilted due to these electromagnetic field 
fluctuations! 

Early time fluctuations in the position of partici- 
pant baryons were discussed in the past as a source 
of the impact parameter fluctuation. Its influence 
survives till the freeze-out resulting in a consid- 
erable difference between the theoretical reaction 
plane and the measured event plane. This effect 
leads to an increase in the magnitude of the elliptic 
flow and generates nonvanishing odd flow harmon- 
ics. 

Wc have found out that within the PHSD model the 
retarded strong electromagnetic field - created dur- 
ing nucleus-nucleus collisions - turns out to be not 
so important as has been expected before. Similarly 
to the HSD results in Ref. Q, the electromagnetic 
field has almost no influence on observables. The 
reason is not a shortness of the interaction time, 
when the electromagnetic field is maximal, but the 
compensation of the mutual action of transverse 
electric and magnetic components. This compen- 
sation effect might be important, for example, if 
an additional induced electric field (as a source of 
the CME) is available in the system since this field 
will not be entangled due to other electromagnetic 
sources. 

Another important point emerging from the com- 
pensation effect of electric and magnetic forces is 
worth mentioning: A significance of an external 
magnetic field in astrophysics is largely accepted. 
There are many studies where various effects of ex- 
ternal magnetic fields are discussed in the appli- 
cation to astrophysics (e.g., see the Introduction 
in Ref. 0] and references in (64|. It is correct in 
this particular problems, however, in many cases it 
is concluded by a statement like "the same effect 
should be observed in high-energy heavy-ion colli- 
sions" which does not hold true due the compensa- 
tion effect as demonstrated in the present work. 

In the intermediate stage of the heavy-ion collision 
the statistical fluctuations of charged quasiparticlcs 
in momentum space can generate charge dipoles 
or even charge quadrupoles. However, the magni- 
tudes Q c i and Q C 2 are small; their orientation is 
distributed almost uniformly and the direction of 
the main axis is changed from event to event. The 
influence of the electromagnetic field here is negli- 
gible again. 



• The transverse momentum conservation - proposed 
as an alternative mechanism for an explanation of 
the observed azimuthal asymmetry - shows a corre- 
lation of the CME and the elliptic flow. However, 
the effect estimated at ^/snn = 200 GcV is too 
small and insensitive to the charge separation. 

• A possible charge separation of balancing charges 
has been addressed by the charge balance function. 
We note that the PHSD model fails to describe the 
focusing effect of the balance function for central 
Au+Au collisions. Certainly, further investigations 
of this problem are needed, both in theory and ex- 
periment especially at lower energies. 

The PHSD approach naturally takes into account the 
main alternative mechanisms of the CME: the momen- 
tum conservation and local charge conservation as well 
as clusters (mini jets, strings, prehadrons, resonances). 
At the moderate energies ^/JWn = 7.7 and 11.5 GeV the 
PHSD model results are close to the experiment since 
partonic degrees-of-freedom are sublcading. However, at 
higher collision energy the PHSD model fails to repro- 
duce the observed azimuthal asymmetry. In contrast 



with our earlier analysis within the HSD model [27j, the 
PHSD overestimates the measured centrality dependence 
of azimuthal distributions due to an increasing action of 
the repulsive scalar parton potential which generates the 
collective flow harmonics in accordance with experiment. 
This finding suggests that a new source of azimuthal 
anisotropy fluctuation is needed beyond the 'standard' 
interactions incorporated in PHSD. The new source does 
not dominate in out-of-plane direction as could be ex- 
pected for the CME but both in-plane and out-of-planc 
components contribute with a comparable strength. In 
this respect the interpretation of the CME STAR mea- 
surements is still puzzling. 

The present PHSD model is already quite elaborated, 
however, as our analysis has shown, color degrees-of- 
freedom or intimate peculiarities of non-Abelian Yang- 
Mills theory should additionally be taken into consid- 
eration. In particular, this concerns the very early 
stage of the nuclear interaction. In this initial state the 
highly compressed strongly interacting matter is dense 
and though the QCD coupling constant is small, glu- 
onic states have high occupation numbers, i.e. the par- 
tons begin to overlap in phase space which leads to some 
saturated state. Strong color forces might create strong 
chromo-electric and chromo-magnetic fields producing a 
new state, a glasma [ill |43| . |65| as mentioned above in 
context of the discussion in Subsection IIII A[ or forming 
new objects like "string ro pes " described in the frame- 
work of Yang-Mills theory [66| . We are planning to in- 
clude these effects into the PHSD model in near future. 

Another class of strong fields relevant to the chirality 
and confinement of dynamical quarks is the long range 
(or soft) vacuum gluon field configurations. Long range 
vacuum gluon fields can be seen as an origin of nozcro 
gluon condensate and topological susceptibility of QCD 
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vacuum [67|, l68| . Soft fields arise in the consideration of 
the global minima of the QCD effective action[69[ and are 
known to play an important role in hadron phenomenol- 
ogy at zero temperature (70j . The nonzero gluon con- 
densate survives at high temperature as demonstrated by 
QCD lattice calculations [7l| . Interplay of strong electro- 
magnetic and vacuum long-range gluon fields can lead to 
the qualitatively new effects in high energy heavy ion col- 
lisions [72l |. However these effects are beyond the scope 
of this paper. 
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